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Justification  is  provided  for  modelling  the  arrival  and 
subsequent  disposition  of  distress  cases  at  a  Coast  Guard 
Search  and  Rescue  Station  as  an  M/E./C  queue.   Numerical 
tables  and  a  simple  computational  method  are  presented  which 
provide  convenient  and  accurate  steady  state  solutions  to 
this  system  assuming  first  in  first  out  queue  discipline. 
Specific  measures  of  effectiveness  computed  include  the 
delay  probability,  average  queue  length,  average  number  in 
system,  average  waiting  time,  average  time  in  system,  and 
the  distribution  of  waiting  time.   Additionally,  the  conse- 
quences of  considering  a  preemptive  priority  queue  discipline 
with  two  priority  classes  are  investigated  and  a  method  of 
calculating  the  delay  probability  for  this  priority  scheme 
is  proposed. 
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I.   INTRODUCTION 

Search  and  Rescue  (SAR)  is  one  of  the  United  States 
Coast  Guard's  most  visible  and  well  known  primary  missions. 
The  backbone  of  the  Coast  Guard's  SAR  effort  is  the  network 
of  small  boat  Stations  distributed  at  various  points  along 
our  nation's  shoreline  so  as  to  provide  adequate  coverage 
for  coastal  boating.   The  bulk  of  these  Stations'  SAR  fleets 
consists  of  two  boat  types,  the  forty  one  foot  utility  boat 
(UTB)  and  the  forty  four  foot  motor  surf  boat  (MSB) .   These 
two  hull  types  account  for  virtually  all  the  short  range 
surface  SAR  response  provided  by  the  Coast  Guard. 

While  it  would  be  ideal  in  terms  of  public  safety  to 
maintain  an  unlimited  number  of  SAR  resources  available  for 
immediate  response,  such  a  policy  is  economically  infeasi- 
ble.   This  gives  rise  to  the  obvious  though  difficult  prob- 
lem of  determining  the  optimal  allocation  of  a  scarce 
resource,  in  this  case,  boat  crews.   The  allocation  of  boat 
crews  vice  boats  is  discussed  since  it  is  generally  Coast 
Guard  policy  to  allow  at  least  one  more  boat  than  boat  crew 
per  Station  to  increase  the  likelihood  of  having  both  a  boat 
and  crew  available  for  response  simultaneously.   Thus,  the 
boat  crew  is  the  limiting  resource. 

The  Search  and  Rescue  Division  of  the  Office  of  Opera- 
tions, Coast  Guard  Headquarters,  maintains  ultimate  responsi- 
bility for  this  task.   In  deciding  on  each  Station's  manning 


level,  Division  personnel  are  interested  in  the  effects 
felt  by  local  SAR  systems  caused  by  varying  available  re- 
sources.  Specific  effects  can  be  measured  by  such  measures 
of  effectiveness  (MOE's)  as  the  probability  of  a  delayed 
distress  response  caused  by  occupied  boat  crews,  the  average 
duration  of  the  delay,  and  the  average  number  of  distressed 
craft  having  to  wait.   Obtaining  such  information  would 
surely  assist  planners  in  arriving  at  some  satisfactory 
allocation . 

These  MOE's  suggest  the  use  of  queueing  techniques  to 
analyze  the  operation  of  SAR  resources  at  a  Station.   The 
incidence  of  distress  cases  can  be  considered  as  customer 
arrivals.   Upon  entering  the  system,  cases  may  either  wait 
for  service  or  receive  service  immediately  upon  finding  an 
idle  boat  crew.   The  length  of  time  a  boat  is  deployed  on  a 
given  case  will  comprise  the  service  time.   Information  on 
both  arrival  times  and  service  times  is  contained  in  the  SAR 
data  base  maintained  at  Coast  Guard  Headquarters.   The  data 
base  consists  of  information  recorded  by  each  unit  on  the 
Search  and  Rescue  Assistance  Report  (CG-5151)  [Ref.  1]  each 
time  they  respond  to  a  case. 

This  thesis  is  concerned  with  examining  the  available 
data  in  order  to  determine  an  appropriate  queueing  system 
which  will  best  model  the  situation,  and  then  providing  solu- 
tions for  the  system  in  terms  of  some  MOE's  which  could  be 
of  value  to  the  SAR  planner.   Chapter  II  of  this  thesis 


contains  a  description  of  the  manner  in  which  distress  cases 
arrive  at  a  Station.   Justification  for  considering  only  a 
finite  time  window  in  order  to  identify  a  period  of  homo- 
geneous arrivals  during  the  busiest  SAR  season  is  presented, 
along  with  the  results  of  the  data-fitting  analysis.   A 
discussion  of  service  times  follows  in  Chapter  III,  again 
along  with  the  results  of  efforts  to  find  an  appropriate 
distribution  to  fit  the  data.   The  samples  were  taken  from 
four  environmentally  diverse  Stations  intended  to  represent 
a  reasonable  cross  section  of  Coast  Guard  operating  areas. 
The  information  was  drawn  from  the  fiscal  year  19  80  data  base 
as  it  was  the  most  recent  available  at  the  time  of  commencing 
this  study.   Chapter  IV  provides  a  list  of  desirable  system 
descriptors  of  use  to  a  Coast  Guard  SAR  planner  and  equates 
these  to  some  common  measures  of  effectiveness  often  used 
in  queueing  analysis.   Two  practical  computational  methods 
for  calculating  the  desired  information  and  some  work  con- 
sidering the  introduction  of  a  two  class  priority  scheme  are 
then  presented  in  Chapter  V.   The  final  chapter  contains  a 
brief  summary  of  the  work. 
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II.   DISCUSSION  OF  THE  ARRIVAL  PROCESS 

A.   DEFINITION  OF  AN  ARRIVAL 

The  definition  of  an  "arrival"  warrants  some  discussion. 
From  the  time  a  vessel  perceives  its  distress  to  the  moment 
a  Coast  Guard  vessel  arrives  alongside,  several  communica- 
tions have  transpired,  both  between  the  distressed  vessel 
and  the  Coast  Guard,  and  within  the  Coast  Guard's  SAR  struc- 
ture.  Exactly  when  during  this  time  a  case  is  said  to 
"arrive"  needs  to  be  decided. 

A  distressed  craft  will  communicate  its  situation  by 
any  means  available.   This  could  be  by  a  mayday  message 
broadcasted  over  a  radio  if  so  equipped,  by  one  of  many 
visual  distress  signals  such  as  lighting  a  fire  on  deck  or 
sending  up  a  red  flare,  or  maybe  simply  by  shouting  to 
another  vessel  in  close  proximity.   This  initial  call  for 
help  may  or  may  not  fall  directly  on  Coast  Guard  ears,  but, 
in  most  cases,  the  Coast  Guard's  first  notification  will 
arrive  at  a  Radio  Station.   While  not  normally  equipped  with 
life  saving  resources,  they  do  act  as  a  communications 
facility  and  can  effectively  relay  relevant  information  to 
the  cognizant  Rescue  Coordination  Center  (RCC).   The  RCC 
acts  as  a  command  center.   Here  the  decision  is  made  as  to 
the  appropriate  resource  to  dedicate  to  the  mission.   If  it 
is  felt  the  best  response  could  be  provided  by  a  Station's 
OTB  or  MSB,  that  Station  is  contacted,  either  by  radio-teletype 
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or  telephone,  and  ordered  to  dispatch  on  the  case.   The 
Station  then  assigns  a  boat  and  a  crew  and  gets  them  under- 
way towards  the  vessel  in  distress. 

Since  only  those  data  recorded  on  the  SAR  Assistance 
Report  [Ref.  1]  are  available,  three  entries  emerge  as 
candidates:   day/time  CG  notified,  day/time  unit  notified, 
and  day/time  underway.   Of  the  three,  the  second,  day/time 
unit  notified,  appears  to  be  the  most  appropriate.   The 
purpose  of  this  work,  is  to  gain  insight  into  the  conse- 
quences of  varying  resources  at  the  Station,  other  factors 
remaining  constant.   It  therefore  would  not  be  appropriate 
to  consider  a  case  as  having  arrived  prior  to  the  Station 
being  made  aware  of  its  need  for  response,  eliminating  the 
first  candidate.   As  for  the  third,  day/time  underway  would 
exclude  time  spent  in  the  queue  awaiting  an  idle  resource, 
while  it  is  an  aim  of  this  thesis  to  use  that  quantity  as 
an  MOE.   This  leaves  day/time  unit  notified  as  the  logical 
time  of  arrival. 

B.   DATA  ANALYSIS 

1.   Homogeneous  Time  Window 

As  suggested  in  the  Introduction,  Station  SAR  demand 
is  far  from  homogeneous,  due  largely  to  the  fact  that  the 
"customers"  serviced  are  primarily  members  of  the  recreational 
boating  population.   Experience  and  the  data  clearly  show 
SAR  case  load  and  boating  density  are  lightest  during  the 
colder,  more  severe  winter  months  and  during  the  normal  five 
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day  work  week.   Rather  than  consider  these  non-homogeneous 
arrivals  in  the  queueing  analysis,  an  attempt  was  made  to 
find  a  finite  time  window  during  the  busiest  hours  of  the 
busiest  days  of  the  busiest  month  in  which  arrivals  could 
be  considered  to  be  occurring  at  a  constant  rate. 
2 .   Statistical  Procedure 

Graphical  displays  of  weekend  SAR  data  suggested  the 
period  0900  to  1800  on  Saturday  and  Sunday  as  a  candidate 
for  a  time  interval  over  which  cases  arrive  at  a  constant 
rate.   Times  from  the  beginning  of  the  time  window  to  the 
first  case  arrival  and  from  the  last  arrival  to  the  end  of 
the  window  were  not  considered.   Efforts  to  fit  these  inter- 
arrival  times,  drawn  from  four  environmentally  diverse 
Stations  (New  Haven,  Ct.,  Rockaway,  N.Y.,  Cape  May,  N.J., 
and  Coos  Bay,  Ore.) ,  to  the  Exponential  distribution  proved 
successful.   July  data  was  chosen  to  represent  the  season 
of  busiest  SAR  activity.   The  decision  to  accept  the  Exponen- 
tial distribution  was  based  on  Chi-square  goodness  of  fit 
tests  at  the  95%  level,  using  the  first  moment  estimator, 
the  reciprocal  of  the  sample  mean  interarrival  time,  as  the 
estimate  of  the  arrival  rate  A.   Table  I  contains  a  summary 
of  the  results  of  the  tests  on  each  Station's  interarrival 
data.   The  computed  Chi-square  statistics  are  all  well  below 
the  corresponding  critical  values,  providing  no  evidence  to 
reject  the  assumption  that  the  data  came  from  an  Exponential 
population.   The  acceptance,  therefore,  of  the  Exponential 
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TABLE  I 
Exponential  Chi-Square  Test  Results 


Station 

X    (hrs.)"1 

2 

^computed 

2 

x .  9  5 , n- 1 

Cape  May 

1.797 

3.30 

12.59 

Coos  Bay 

1.022 

2.37 

12.59 

New  Haven 

1.018 

0.35 

9.49 

Rock aw ay 

1.340 

3.22 

14.07 

distribution  implies  that  the  Poisson  Process  is  an  appro- 
priate descriptor  of  the  arrival  phenomenon. 
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III.   DISCUSSION  OF  SERVICE  TIMES 

A.   DEFINITION  OF  SERVICE  TIME 

As  with  arrivals,  exactly  what  is  to  be  considered 
service  time  needs  to  be  defined.   Unlike  a  bank  or  barber 
shop,  customers  do  not  physically  arrive  at  a  well-defined 
facility  and  commence  being  administered  to,  face  to  face, 
by  the  first  available  server,  nor  do  customers  "queue  up" 
in  a  waiting  line  when  awaiting  service.   In  most  instances, 
a  typical  SAR  case  occurs  some  distance  offshore.   The  crew 
must  get  a  normally  moored  vessel  underway,  transit  the 
distance  from  the  Station  to  the  distressed  craft,  render 
necessary  assistance,  and,  in  more  cases  than  not,  tow  the 
disabled  vessel  to  safety.   Since  all  elements  of  the  afore- 
mentioned procedure  are  part  of  servicing  a  case,  all  should 
contribute  to  service  time. 

A  considerable  problem  arises  due  to  the  nature  of  the 
available  information.   Data  recorded  include  day/time  unit 
notified  and  day/time  underway  as  previously  discussed,  as 
well  as  sortie  time.   This  latter  term  is  defined  in  the 
Coast  Guard  publication  Search  and  Rescue  Data  System  Manual 
(COMDTINST  M5230.10)  [Ref.  2]  and  basically  involves  the 
underway  time  of  the  rescue  vessel.   More  specifically, 
sortie  time  is  either  the  time  from  departure  until  the  time 
of  return  to  homeport ,  or,  if  diverted  from  another  mission, 
the  time  of  diversion  to  the  time  of  return  to  the  original 
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mission.   The  problem  arises  in  considering  the  elapsed  time 
from  notification  to  time  underway.   Given  an  available 
resource  (boat) ,  this  time  typically  involves  either  a  short 
walk  from  an  adjacent  shore  facility  to  the  boat  or  often  a 
radio  call  directly  to  the  boat  if  the  crew  is  already 
aboard.   However,  if  a  case  were  queued,  the  waiting  time 
would  accumulate  here  with  no  accompanying  indication  that  a 
resource  was  unavailable.   In  other  words,  there  is  no  way  to 
isolate  that  portion  of  the  time  from  notification  to  the 
time  underway  due  to  normal  startup  time,  which  of  course 
should  be  included  in  service  time,  from  that  time  spent 
waiting  for  the  availability  of  a  resource.   Examination  of 
the  raw  data,  however,  revealed  that  the  time  from  unit 
notification  to  time  underway,  which  would  include  any 
queueing  time,  was  no  more  than  ten  minutes  90%  of  the  time. 
In  further  analysis,  therefore,  sortie  time  will  be  used  as 
a  surrogate  of  service  time. 

B.   DATA  ANALYSIS 

1 .   Parameter  Estimation 

Sortie  times  for  the  four  Stations  were  then  dis- 
played in  stemleaf  plots  to  aid  in  discovering  an  appropriate 
distribution  to  fit  the  data.  Only  those  service  periods 
begun  during  the  homogeneous  time  window  were  considered. 
Examination  of  the  displays  consistently  suggested  the  use 
of  the  Gamma  or  Erlang  distributions,  of  which  the  latter 
was  chosen  due  to  the  availability  of  many  queueing  results 
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involving  its  use.   Due  to  the  presence  of  two  parameters  in 
the  Erlang  distribution,  however,  three  methods  of  parameter 
estimation  were  investigated.   Letting  Y  be  a  random  variable 
representing  service  time,  the  appropriate  PDF  is 


.1,K  1 

f(y)   =   (kt1}  ,  Y    e  "   ,    y  >  0  (3.1) 


2 
with  mean  k3  and  variance  k3  . 

The  first  alternative  involved  using  the  method  of 
moments,  setting  the  sample  moments  equal  to  the  correspond- 
ing distributional  moments; 


k8      and      S2   =   kg2. 


This  yielded  the  parameter  estimates: 


k   =   — j  and      3   =   — , 

S^  Y 


Since  k  must  take  on  an  integer  value,  k  must  be  rounded 
to  the  nearest  integer. 

The  second  method  varies  from  the  first  only  in  its 
estimate  for  3  and  accounts  for  the  rounding  error  created 
when  rounding  k  to  the  nearest  integer.   The  estimate  for 
3  is  therefore  recalculated  using  the  integer  valued  estimate 
of  k  and  the  relationship 
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mean   =  kg 


yielding 


Y 

3   =  r  • 


The  third  method  investigated  is  presented  by  Read 
[Ref .  3]  and  is  derived  through  the  Gamma  distribution,  then 
making  the  appropriate  conversions  to  the  Erlang  distribu- 
tion.  His  estimates  for  the  Gamma  parameters 


Yz  ,       ~      -   1 

a   =   and      S   =   Y  -  — 


Yz-1 


z 


where 


z   = 


1   n   l 

n   £   Y. 

1=1   l 


are  said  to  be  asymptotically  more  efficient  than  the  simple 
method  of  moments  estimators  for  k  >  3. 
2 .   Statistical  Procedure 

The  Erlang  parameter  estimates  for  each  of  the  four 
Station  samples  were  calculated  using  each  of  the  aforemen- 
tioned methods.   Using  the  estimates  provided,  Chi-square 
goodness  of  fit  tests  were  conducted.   The  second  method 
consistently  produced  the  smallest  Chi-square  statistics, 
well  within  a  95%  acceptance  region,  and  thus  was  chosen  to 
provide  future  parametric  estimates.   The  results  are 
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displayed  in  Table  II.   It  should  be  noted  that  the  perfor- 
mance of  the  third  method  in  terms  of  the  Chi-square  tests 
was  not  surprising  due  to  the  small  values  of  k  generated. 
Perhaps  in  future  application,  if  samples  produce  larger 
estimates  of  k,  the  use  of  this  method  would  prove  to  be  a 
more  competitive  alternative  and  should  therefore  not  be 
completely  abandoned. 

TABLE  II 
Erlang  Chi-Square  Test  Results 
Method  #1        Method  #2        Method  #3 


Station 

k 

e(hrs.) 

2 

y 

Acomp 

k 

B(hrs.) 

2 
\x>mp 

k 

3(hrs.) 

2 

Aoomp 

2 

x.95,n-l 

Cape  May 

2 

.703 

10.16 

2 

.778 

8.24 

2 

.668 

21.01 

12.59 

Coos  Bay 

1 

1.169 

4.27 

1 

1.051 

2.18 

2 

.573 

11.24 

9.49 

New  Haven 

2 

.810 

11.84 

2 

.689 

6.70 

3 

.534 

11.16 

9.49 

Rockaway 

2 

.620 

7.39 

2 

.577 

6.14 

2 

.611 

5.39 

11.07 

19 


IV.   MEASURES  OF  EFFECTIVENESS 

The  results  from  the  two  previous  chapters  support  the 
contention  that  the  M/E, /c  queue  may  be  of  value  in  des- 
cribing the  operation  and  allocation  of  Search  and  Rescue 
resources  at  a  Coast  Guard  Station.   It  remains  to  determine 
a  method  of  extracting  needed  measures  of  effectiveness  to 
aid  in  evaluating  the  relative  merits  of  competing  resource 
allocations. 

In  administering  the  Search  and  Rescue  program,  Coast 
Guard  planners  must  decide  upon,  as  previously  discussed,  the 
number  of  ready  boat  crews  to  allocate  to  each  Station. 
Their  problem  is  that  of  supplying  the  customers,  in  this 
case  the  recreational  boating  public,  with  some  satisfactory 
level  of  SAR  capability  while  attempting  to  minimize  system 
costs.   Certain  relevant  measures  of  system  performance 
obtainable  through  queueing  techniques  could  prove  of  great 
utility  in  this  endeavor.   These  measures  are  described 
below . 

A.   PROBABILITY  OF  A  DISTRESS  CASE  WAITING  FOR  SAR  SERVICE 

An  important  yardstick  with  which  to  measure  the  adequacy 
of  a  given  resource  level  is  the  estimate  of  the  probability 
that  a  vessel  in  distress  will  have  to  wait  for  an  idle  res- 
cue craft.   Due  to  the  emergency  nature  of  the  required 
response,  it  is  of  paramount  importance  to  maintain  a  con- 
sistently high  readiness  posture.   Unlike  a  barber  shop  or 
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gas  station,  the  consequences  of  waiting  for  service  may  be 
grave  and  even  potentially  fatal.   A  method  through  which 
one  can  assess,  for  varying  boat  crew  allocations,  the  likeli- 
hood of  failure  to  meet  the  SAR  demand  at  once  during  a 
peak  period  would  be  of  great  value.   Using  queueing  analy- 
sis, the  quantity  termed  delay  probability  provides  this 
result  and  will  be  calculated  in  the  next  chapter. 

B.  AVERAGE  NUMBER  OF  DISTRESS  CASES  WAITING  FOR  SERVICE 
AND  IN  THE  SYSTEM 

The  analyst  might  choose  to  weigh  the  relative  merits  of 
two  or  more  proposed  allocation  schemes  based  on  the  number 
of  distressed  vessels,  on  the  average,  having  to  wait  for  an 
idle  Coast  Guard  crew,  or  similarly,  on  the  average  number 
of  active  cases,  both  those  being  assisted  and  those  waiting 
for  service,  at  any  time.   Clearly,  the  ordinal  ranking  of 
fewer  is  better  would  be  appropriate.   The  queueing  expres- 
sions L   (Average  queue  length)  and  L  (average  number  in 

Si 

system)  provide  methods  of  calculating  these  averages. 

C.  AVERAGE  TIME  SPENT  WAITING  BY  A  DISTRESS  CASE  FOR  AN 
IDLE  CREW  AND  AVERAGE  TIME  IN  SYSTEM 

The  planner  might  wish  to  provide  boat  crews  sufficient 

to  keep  the  average  duration  of  a  distressed  case's  delay 

below  some  maximum  value.   The  total  time,  on  the  average,  a 

case  remains  active,  including  both  waiting  and  service  time, 

might  provide  another  MOE  of  similar  value.   These  quantities, 

designated  W   (average  waiting  time)  and  W  (average  time  in 

si 

system)  are  readily  obtainable  through  queueing  techniques. 
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D.   DISTRIBUTION  OF  WAITING  TIMES  OF  A  DISTRESS  CASE 

Often  a  SAR  system  manager  would  be  interested  in  the 
likelihood  that  a  distressed  vessel's  delay  exceeds  some 
fixed  length  of  time.   This  type  of  analysis  lends  itself  to 
the  establishment  of  some  threshold  duration  or  maximum 
acceptable  delay  with  very  high  probability.   Queueing 
analysis  could  allow  the  planner  to  determine,  for  instance, 
the  number  of  boat  crews  required  to  ensure  that  the  likeli- 
hood a  vessel  must  wait  more  than  a  half  hour  for  an  idle 
rescue  boat  does  not  exceed  ten  percent.   A  knowledge  of  the 
distribution  of  waiting  times  would  therefore  serve  as 
excellent  input  to  the  decision  making  process. 
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V.   COMPUTATIONAL  METHODS 

This  chapter  presents  two  reasonably  simple  methods  for 
obtaining  the  aforementioned  MOE ' s  given  the  appropriate- 
ness of  the  M/E,/c  queue  as  a  modelling  vehicle  for  the 
Station  SAR  operation,  and  then  extends  the  results  to 
consider  a  priority  scheme  other  than  FIFO. 

The  work  thus  far  has  concentrated  on  providing  neces- 
sary background  on  the  Coast  Guard  SAR  system,  finding 
characterizations  of  the  arrival  process  and  service  times, 
and  selecting  MOE ' s .   Prior  to  presenting  the  two  computa- 
tional methods  previously  discussed,  a  more  detailed  descrip- 
tion of  the  M/E, /c  queueing  system  is  in  order. 

A.   DESCRIPTION  OF  THE  M/Ek/C  QUEUE 

The  M/E,  /c  queue  is  characterized  by  a  homogeneous 
Poisson  arrival  stream  with  rate  A  whose  interarrival  times 
have  PDF 


f(x)   =  X   e"Xx,    x  >  0 


and  Erlang  service  times  with  PDF 


(-)k        -- 

f(y)  =  jy^jr  yk_1  e  sY,   y  i  o 
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The  FIFO  (first  in,  first  out)  priority  discipline  is 
assumed  and  the  following  variable  and  parameter  definitions 
apply. 

c   =   total  #  of  servers 

A    =   arrival  rate   =   reciprocal  of  mean  inter- 
arrival  time 

k   =   Erlang  shape  parameter 

S    =   Erlang  scale  parameter 

U    =   service  rate   =   reciprocal  of  mean  service 
time 

p    =   A/cy   =   system  intensity 

N   =   total  number  in  system 

N    =   total  number  in  queue 

q 

L   =  E[N]   =   mean  number  in  system 
L   =  E [N  ]   =  mean  number  in  queue 

q      q 

T   =   total  time  in  system 
T   =   total  time  in  queue 

q 

W   =   E[T]   =   mean  time  in  system 
W    =   E [T  ]   =   mean  time  in  aueue 

q      q 

B.   TABULAR  METHOD 
1 .   Background 

The  theory  surrounding  the  M/E,  /c  queue  has  been 
sufficiently  developed  to  a  point  where  several  numerical 
methods  have  been  shown  to  produce  adequate  steady  state 
results.   Due  to  the  complexity  of  the  computational  require- 
ments, however,  direct  utilization  by  a  layman  in  a  decision 
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making  capacity  (in  this  case,  a  Lieutenant  in  the  Search 
and  Rescue  Branch)  would  be  impractical,  if  not  impossible. 
Hillier  and  Lo  [Ref.  4],  in  an  effort  to  make  complex  results 
readily  obtainable  by  would-be  users,  have  prepared,  for  a 
finite  though  practical  subset  of  possible  values  of  system 
parameters  k,  c,  and  p,  a  series  of  tables  presenting  many 
desirable  measures  of  steady  state  system  performance.   They 
cite  recent  computational  treatments  by  Mayhugh  and  McCormick 
[Ref.  6],  Mayhugh  [Ref.  7],  and  Heffer  [Ref.  8],  and  the 
comprehensive  analysis  of  the  E  /E,/c  system  by  Yu  [Ref.  9] 
as  the  basis  for  their  tables.   Of  specific  interest  for 
this  application  is  their  Table  I,  partially  reproduced 
herein  as  Appendix  A,  which  provides  the  desired  measures  of 
effectiveness  for  the  following  parametric  combinations. 

k  =  2  ;   c   =   2  ,  3  ,  .  .  .  ,  10 

k  =  3;   c   =   2,3,4,5 

k  =  4 ;   c   =   2,3 

k  =  5,6,7,8;   k   =   2 

g  =  .10, .20, . . . ,  .50,  .55, . . . , .95, .98, .99 

2 .   Use  of  Tables 

Use  of  the  tables  requires  entering  arguments  c,  k, 
and  p.   The  parameter  c,  number  of  boat  crews,  will  most 
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likely  be  the  design  variable  whose  determination  is  the 
object  of  analysis.   Estimates  for  the  parameters  k  and 
p  will  be  provided  by  the  SAR  data  from  the  Station  under 
study.   Recommended  parametric  estimates,  as  previously 
discussed,  are 

X      =   1/X 

—2   2 

k   =   Y  /S    then  rounding  to  nearest  integer 

3   =   Y/k 

/\  A  A 

U  =  IAS 

p   =   A/cy 

The  reader  is  referred  to  Chapters  II  and  III  for  discussions 
concerning  the  applicability  of  the  SAR  data  base  and  the 
criteria  for  the  determination  of  the  homogeneous  time  window 
to  satisfy  the  assumptions  of  the  M/E^/c  queue. 

Since  results  are  listed  for  discrete  values  of  the 
system  intensity  p,  interpolation  will  undoubtedly  be  neces- 
sary.  Since  no  interpolation  scheme  is  suggested  in  Ref.  4, 
and  since  visual  examination  of  the  tabular  entries  suggests 
it  would  not  be  unreasonable,  the  linear  technique  is  assumed, 

The  aforementioned  measures  of  effectiveness  are 
obtainable,  directly  or  indirectly,  from  the  tables.   The 
probability  of  requiring  a  distress  case  to  wait  for  response 
due  to  preoccupied  rescue  boats,  otherwise  called  delay 
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probability,  is  listed  separately  at  the  end  of  each  section 
of  parametric  combinations.   Though  interpolation  for  the 
exact  value  of  system  intensity  p  is  required,  this  impor- 
tant result  is  otherwise  available  at  the  analyst's  finger- 
tips.  Likewise,  average  number  of  cases  waiting  (L  )  and 
average  number  of  cases  in  the  system  (L)  are  similarly 
accessible.   Calculation  of  average  time  spent  waiting  for 
service  (W  )  and  average  time  in  the  system  (W)  is  made 

q 

possible  through  the  well  known  Little's  [Ref.  10]  relation- 
ships 


L 

W   =   -t3        and        W   =   -  . 

q      A  X 


In  addition,  though  not  directly  tabulated,  Hillier 
suggests  an  approximation  method  for  estimating  the  distri- 
bution of  waiting  time  T  .   The  derivation  is  based  upon  the 
observation  that  an  Erlang  random  variable  with  mean  1/n  may 
be  expressed  as  the  sum  of  k  Exponential  phases,  each  with 
mean  1/ky,  and  that,  when  there  are  no  idle  servers,  comple- 
tions of  service  phases  occur  according  to  a  Poisson  Process 
with  rate  cky .   This  leads,  as  discussed  in  Ref.  4,  to  the 
following  probability  distribution  function. 


P  [T  >  t]   =       P  [N=n]P[D  <_k(n-c+l) -1]  ,   t  >_  0 
"  n=c 


where  D  is  a  Poisson  random  variable  with  mean  ckut.   Since 
the  first  term  is  tabulated  by  Hillier  and  Lo  in  Ref.  4,  and 
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cumulative  Poisson  probabilities  are  readily  obtainable, 
either  in  tabular  form  or  through  direct  computation,  the 
entire  expression  is  easily  computable. 
3.   Examples 

By  manner  of  illustration,  the  following  two  tables 
demonstrate  the  effects  on  system  behavior  resulting  from 
varying  the  number  of  servers  (boat  crews)  in  an  otherwise 
stable  queue.   Table  III  contains  results  stemming  from  the 
author's  arbitrary  and  hypothetical  selection  of  system 
parameters  A  and  y,  while  Table  IV  displays  results  produced 
from  actual  SAR  data  drawn  from  Coast  Guard  Station  Rockaway, 
N.Y. 

TABLE  III 
Hypothetical  Case 

A  =  1.5/hr    k  =  3    \i   -  1.0/hr 


c   =   2 

c   =    3 

c    =    5 

MOE 

p    =    .75 

p    =    .50 

p    =    .30 

P[N  =  0] 

.139 

.208 

.223 

P[T     >  0] 

q 

.6  39 

.232 

.020 

L 

2.803 

1.667 

1.507 

Lq 

1.30  3 

.167 

.007 

W    (hrs.) 

1.869 

1.111 

1.005 

W      (hrs.) 

q 

.869 

.111 

.005 

P  [T     >   . 5 

q 

hrs.  ] 

.513 

.105 

.004 
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TABLE  IV 
Station  Rockaway 


X    =  1.340/hr 


k  =  2 


B  =  .577  hrs 


y  =  .86  7/hr 


c   =    2 

c   =    3 

c 

=    4 

MOE 

p      =    .77 

p    =    .52 

P 

=    .39 

P[N  =  0] 

.128 

.195 

.209 

P  [T     >  0] 

q 

.668 

.256 

.084 

L 

3.275 

1.782 

1.605 

Lq 

1.735 

.222 

.045 

W    (hrs.) 

2.444 

1.330 

1.197 

W      (hrs.) 

q 

1.295 

.165 

.033 

P[T     >   .5 

q 

hrs.  ] 

.561 

.088 

.008 

The  above  tables  illustrate  quite  graphically  the 
effects  of  varying  the  number  of  servers  or  boat  crews  and 
could  justifiably  provide  the  SAR  planner  with  valuable 
information  upon  which  to  base  a  decision.   Given  the  exis- 
tence of  some  design  criterion  specifying,  for  example,  a 
maximum  acceptable  delay  probability,  this  method  of  analysis 
would  clearly  indicate  the  minimum  number  of  boats  required 
to  achieve  it.   For  example,  if  Coast  Guard  policy  required 
that  the  chance  of  a  delayed  response  not  exceed,  say,  .30, 
Table  IV  clearly  indicates  at  least  three  boat  crews  would 
be  required  at  Station  Rockaway  in  order  to  meet  that  require- 
ment during  the  peak  load  periods.   Similarly,  a  system  re- 
quirement that  the  average  amount  of  time  spent  waiting  for 
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service  be  less  than  some  amount,  say  one  half  hour,  could 
very  simply  be  satisfied  through  ensuring  that  the  required 
number  of  boat  crews  are  assigned. 

C.   COMPUTATIONAL  METHOD 

1 .  Background 

In  his  Master's  Thesis,  Celayir  [Ref.  5]  developed  a 
relatively  simple  analytic  solution  to  the  M/E,/c  queue, 
which  has  the  advantage  of  supplying  rapid,  computable 
results  without  the  necessity  for  extensive  tables  or 
sophisticated  computational  capability.   Additionally,  the 
method  is  not  constrained  to  the  discrete  set  of  system 
parameters  k  and  c  covered  in  Lo  [Ref.  4],  and  thus  provides 
a  valuable  alternative  when  existing  SAR  data  produce  param- 
eter estimates  beyond  the  values  mentioned  in  the  previous 
section.   Further,  Celayir' s  results  were  made  plausible 
through  comparisons  made  with  three  other  sources,  specifically, 
simulation  output,  Hillier  and  Lo ' s  [Ref.  4]  tables,  and 
observed  similarities  with  the  tractable  M/M/c  queue.   The 
reader  is  referred  to  Ref.  5  for  a  more  in-depth  discussion 
of  the  derivation. 

2 .  Computational  Procedure 

Like  the  previous  method,  system  parameters  A,  k,  3, 
U,  and  p  are  required  for  the  calculations  and  should  be 
estimated  from  the  SAR  data  base  as  discussed  in  Chapters 
II  and  III.   One  additional  parameter,  offered  load,  is 
mentioned  by  Celayir  and  is  defined  as  the  ratio  of  the 


30 


arrival  rate  A  to  the  service  rate  y  and  can  be  expressed 
as  the  product  of  the  number  of  crews  c  and  the  system 
intensity  p  yielding 


X 
a      =   —  =   c  p  . 

y 


It  was  observed  by  Celayir  that,  for  equal  values 

of  system  intensity  p,  delay  probabilities  (P[T  >0])  for 

si 

the  M/E, /c  queue  were  virtually  identical  to  those  of  the 
M/M/c  queue.   Thus,  the  delay  probability  for  the  latter 
system  provides  an  excellent  approximation  for  that  of  the 
former  and  is  computable  [Ref.  11]  as 


P'V°)   "   cMc^-'I)  PO  'S"1' 


where 


1  ,  X .  n    1  .  A .  c  /  cu  m-1 


po   -    LL  ct(?   ♦  ar(?  <35£x*  > 


n=0 


This  expression,  though  in  a  different  form,  is  algebraically 
equivalent  to  that  presented  in  Ref.  5.   Should  this  calcu- 
lation prove  too  cumbersome,  as  may  be  the  case  for  c  >  4 
(although,  in  this  application,  a  Station  is  rarely  billeted 
for  more  than  three  boat  crews) ,  Appendix  B  contains  two 
charts  which  provide  delay  probabilities  for  the  M/M/c 
queue  given  the  offered  load  a  as  the  entering  argument. 
Either  of  these  two  methods  provides  a  convenient  and  accurate 


31 


approximation,  without  the  need  for  interpolation,  for  the 
probability  a  distress  case  will  need  to  wait  for  an  idle 
boat  crew. 

An  approximation  for  the  average  number  of  vessels 
waiting  for  service,  L  ,  is  derived  in  Ref.  5  and  is  based 
upon  extensions  of  the  theoretical  work  in  Lo  [Ref.  4].   The 
formula  stems  from  Hillier  and  Lo's  approximation 


Lq   =  i(l  +  g)  (1  ♦  |)  Lq 


where  L    is  the  average  queue  length  for  the  corresponding 

ql 
M/M/c  queue  and  the  function  g  is  a  multiplicative  correction 

factor  developed  by  Lo  [Ref.  4]  to  be  applied  in  the  compu- 
tation of  Lp  when  extrapolating  the  tabular  results  for 
large  values  of  the  parameters  k  or  c.   But,  for  smaller 
values  of  system  parameters  k  and  c,  g  tends  to  zero.   Noting 
this  and  the  fact  that  Lq,  may  be  simply  expressed  as  a 
function  of  the  previously  derived  delay  probability,  that  is 

.  (A/y)C\y 

ql      (c-1) ! (cp-X)z 

=  rr^-r    '    P[T  >  0]  , 

cu-A       q 

the  average  queue  length  may  then  be  written  as 


L    -   i(l  +h   (  — K—)   *P[T  >  0]  . 
q      2     K   cy-A      q 
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It  then  follows,  again  using  Little's  [Ref.  10]  relationship, 
that  the  average  time  spent  waiting  for  an  available  crew 
equals 


L 
W    =  -Z 

q     a 


=    T(1  +  £><cFT>"p*V01 


The  average  number  L  of  cases  in  the  system  at  any  time  and 
the  average  time  W  spent  in  the  system  are  then  expressable 
as 


L   =   L   +  i 

q   y 


and 


w   =   w  +  -  =   t-  . 

q   vi    a 


Thus,  given  system  parameters  A,  k,  y ,  and  c  and  having 
obtained  the  delay  probability  P [T  >0],  approximations  for 
average  number  waiting,  average  time  spent  waiting,  average 
number  in  system,  and  average  time  in  system  are  all  easily 
computable. 

Finally,  Celayir  [Ref.  5]  proposes  a  direct  method 
of  calculating  the  probability  that  the  duration  of  a  dis- 
tressed vessel's  wait  exceeds  a  specified  length  of  time. 
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He  suggests  that  the  cumulative  distribution  function  (CDF) 
of  waiting  times  T   is 

*      q 


FT  (t)   =   1  -  P[T  >  0]  e  bt  ,    t  >  0 

q  q 

where 


h   =   2  c  y  LI.  "  P } 


The  interested  reader  is  referred  to  his  work  for  an  explana- 
tion of  the  derivation.   One  should  note  however,  that  all 
component  parts  of  the  expression  are  obtainable  as  previously 
described.   If  the  SAR  planner,  for  example,  wished  to  know 
the  probability  that  a  potential  case  would  have  to  wait 
longer  than  one  hour  for  an  available  Coast  Guard  response, 
he  would  need  only  to  compute 


P  [T  >  0]  e 

q 


Similarly,  should  he  wish  to  know  the  chances  that  an  arriving 
case  might  wait  longer  than  one  hour  but  less  than  two,  the 
expression 


„  rm    ~  n  ,  -b     - 2b, 
P  [T  >  0  ]  (e    -  e    ) 
4 


would  provide  the  answer.   It  should  be  emphasized  that, 
while  Celayir  [Ref.  5]  showed  that  results  obtained  using 
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this  expression  compared  quite  favorably  with  those  from  his 
simulation  and  from  Hillier  and  Lo's  tables,  considerably 
less  effort  was  required  to  compute  his  formulas. 
3.   Examples 

The  following  two  examples  are  presented  to  provide 
a  comparison  with  Tables  III  and  IV.   All  relevant  system 
parameters  were  held  constant. 

TABLE  V 
Hypothetical  Case 


A 

= 

1/5/hr. 

k    =    3 

c   =    2 

p    =    .75 

y    =    1.0/hr. 

c   =    3 
p    =    .50 

c    =    5 
p    =    .30 

MOE 

a    =    1.50 

a    =    1.50 

a   =    1.50 

P[T     >  0] 

q 

.642 

.238 

.019 

L 

2.784 

1.659 

1.505 

L 

q 

1.284 

.159 

.005 

W    (hrs.) 

1.856 

1.106 

1.003 

W      (hrs.) 

q 

.856 

.106 

.003 

P[T     >  .5 

hrs. 

.] 

.441 

.077 

.001 

35 


TABLE  VI 
Station  Rockaway 
A  =  1.340/hr.     k  =  2     3  =  .577  hrs.     y  =  .867/hr. 


c   =    2 

c   =    3 

c   =    4 

p    =    .77 

p    =    .52 

p    =    .39 

MOE 

a    =    1.54 

a   =    1.56 

a    =    1.56 

P[T     >  0] 

q 

.669 

.254 

.077 

L 

3.252 

1.749 

1.582 

Lq 

1.706 

.203 

.036 

W    (hrs.) 

2.427 

1.305 

1.181 

W      (hrs.) 

q 

1.273 

.151 

.027 

P[T     >  .5 

q 

hrs .  ] 

.513 

.195 

.019 

As  can  be  seen,  Tables  V  and  VI  resulting  from  calcu- 
lations using  Celayir's  method  compare  quite  favorably  with 
the  figures  in  Tables  III  and  IV.   Certainly,  in  the  context 
of  providing  the  SAR  planner  with  a  guide  with  which  to  com- 
pare the  relative  performance  of  a  Search  and  Rescue  Station 
given  various  numbers  of  boat  crews,  both  methods  are  of 
great  utility. 

D.   PRIORITY  CONSIDERATIONS 

The  reader  is  reminded  that  the  queueing  techniques 
applied  thus  far  have  assumed  the  FIFO  priority  discipline. 
This  assumption  may  carry  with  it  far  reaching  ramifications, 
for,  in  practice,  resources  are  assigned  to  distress  cases 
according  to  a  continuous  priority  scale.   Each  arriving 
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case  is  evaluated  and  assessed  as  to  its  severity  relative 
to  other  active  cases.   If,  in  the  judgement  of  the  cogni- 
zant Rescue  Coordination  Center,  the  newly  arrived  distress 
case  provides  a  greater  or  more  imminent  threat  to  life  or 
property  than  other  vessels  currently  being  assisted,  ongoing 
assistance  will  be  terminated  and  that  resource  redirected 
towards  the  more  urgent  situation.   A  possible  scenario  fit- 
ting this  description  might  be  a  Coast  Guard  ULB  diverting 
from  retrieving  a  disabled  pleasure  craft  at  anchor  in  a 
protected  waterway  in  favor  of  rendering  assistance  to  a 
vessel  on  fire  or  in  danger  of  going  aground. 
1.   Background 

Extensive  theoretical  work  concerning  the  inclusion 
of  priority  schemes  other  than  first  come  first  served 
exists,  but  most  is  concentrated  on  the  more  common  M/M/l 
queue,  that  is  single  server  with  Exponential  service. 
Extensions  to  systems  involving  multiple  servers  and  service 
distributions  other  than  Exponential  are  rare  as  well  as 
extremely  complex,  and  clearly  stand  as  an  area  for  further 
research.   The  most  thorough  treatment  encountered  by  the 
author  can  be  found  in  Jaiswal  [Ref.  12],  and  the  reader  is 
referred  there  for  a  proper  initiation  into  the  subject.   By 
way  of  an  overview,  however,  the  following  discussion  is 
provided. 

When  considering  priority  schemes  other  than  FIFO, 
arrivals  are  usually  thought  to  belong  to  one  of  several 
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priority  classes.   Prior  to  the  saturation  of  idle  servers, 
customers  are  served  as  they  arrive.   Once  all  servers  are 
busy,  however,  customers  assume  a  place  in  the  queue  accord- 
ing to  their  priority  class,  the  customers  belonging  to  the 
higher  priority  classes  moving  ahead  of  those  in  lower  classes. 
It  can  be  seen,  then,  that  the  potential  exists  for  customers 
of  a  low  priority  class  to  be  bumped  indefinitely.   The 
nature  of  the  priority  discipline  may  be  further  refined  by 
considering  whether  priorities  are  preemptive  or  non-preemptive 
In  the  latter,  an  arriving  customer  belonging  to  say  the 
highest  priority  class  will  move  to  the  front  of  the  queue 
and  wait  for  the  next  available  server  even  if  customers  of 
a  lower  class  are  currently  in  service.   If,  however,  the 
priority  scheme  is  said  to  be  preemptive,  current  service  on 
a  lower  priority  customer  will  be  terminated  upon  the  arrival 
of  a  customer  belonging  to  a  higher  priority  class.   In  the 
context  of  a  Coast  Guard  SAR  Station,  distress  cases  would 
be  considered  as  having  preemptive  priorities  with  the  capa- 
bility of  interrupting  ongoing  assistance  to  a  case  belonging 
to  a  lower  class. 

2 .   Two  Priority  Class  Case 

It  has  been  suggested  by  personnel  from  the  Search 
and  Rescue  Branch  that  two  priority  classes ,  labeled  serious 
and  routine,  would  suffice.   This  designation  arises  from 
information  regarding  the  severity  of  a  case  recorded,  after 
the  fact,  on  Ref.  1,  and  therefore  contained  in  the  SAR  data 
base.   Thus,  separate  arrival  rates  and  candidate  descriptive 
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distributions  for  both  classes  are  readily  obtainable.   While 
this  type  of  analysis  is  a  gross  approximation  to  the  con- 
tinuous set  of  actual  priority  classes,  it  would  probably 
provide  an  improvement  over  considering  the  FIFO  discipline 
alone. 

While  the  complexities  introduced  by  priority  con- 
siderations create  a  considerable  barrier  in  obtaining 
exact  solutions  to  the  M/E,/c  queue,  one  important  measure 
of  effectiveness,  the  delay  probability,  is  somewhat  tracta- 
ble.  If  one  were  to  restrict  the  number  of  priority  classes 
to  two,  for  instance,  estimates  for  the  delay  probability 
for  both  are  derived  for  the  cases  c  =  1  and  c  =  2. 

Considering  first  the  case  of  the  higher  priority  or 
serious  class  of  distress  cases,  the  direct  application  of 
techniques  discussed  in  Chapters  II,  III,  and  V  is  in  order. 
This  stems  from  the  observation  that,  given  the  preemptive 
priority  discipline,  a  server  occupied  with  a  routine  case 
is  identical,  from  the  point  of  view  of  an  arriving  serious 
case,  to  an  idle  server.   Thus,  the  analyst  could  determine 
the  delay  probability  of  a  serious  case  by  considering  only 
the  arrivals  and  service  times  of  cases  belonging  to  that 
class  from  the  initial  data  base,  and  then  proceeding  as 
before,  while  ignoring  all  routine  cases.   This  would  involve 
fitting  this  abbreviated  data  to  arrival  and  service  distri- 
butions, and  then  calculating  the  desired  quantities  from 
the  resultant  queueing  system.   Although  unable  to  confirm 
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it  (available  data  contains  no  more  than  six  serious  arrivals 
per  Station) ,  the  author  has  seen  no  evidence  to  suggest 
that  the  M/E, /c  system  would  not  again  be  appropriate. 

The  question  of  the  delay  probability  for  the  non- 
serious  or  routine  case  presents  a  somewhat  more  complicated 
problem  arising  from  the  fact  that  a  routine  case  might  be 
delayed  even  after  its  service  has  begun.   Specifically,  if 
all  servers  are  busy,  an  arriving  serious  case  could  preempt 
service  on  a  routine  case  with  the  latter  assuming  a  place 
in  the  queue  at  the  front  of  other  cases  belonging  to  the 
same  class.   The  derivations  that  follow  require  the  intro- 
duction of  some  new  notation. 

D   =   total  delay  time  for  a  routine  case; 

D,   =   initial  delay  encountered  by  a  routine  case 
upon  arrival  before  its  service  begins; 

D~   =   subsequent  delay  due  to  one  or  more  preemp- 
tions from  service  by  a  newly  arriving  serious 
case  ; 

N   =   number  in  system  at  time  of  routine  case's 
arrival; 

NR   =   number  of  routine  cases  in  system  at  time  of 
new  routine  case's  arrival; 

A~   =   number  of  new  serious  cases  arriving  during 
service  time  of  a  routine  case. 

a.   Single  Server  System 

The  case  where  there  is  only  one  server  available 
is  considered  first.   Noting  that 


D   =   Dx  +  D2 
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for  a  routine  case,  the  chance  of  delay  can  be  expressed  as 


P[D>  0]   =   1  -  P[D=  0] 


1  -  P[D,  +  D0  =  0] 


1  -  P[D1  =  0,D2  =  0 


=  l  -  p[d2  =  o|d1  =  o]  P[D1  =  0] 


Since  the  initial  delay  D,  is  exactly  the  delay  that  may 
occur  in  a  FIFO  (one  priority)  system,  P[D,  =0]  is  simply 
one  minus  the  delay  probability,  P [T  >  0],  as  calculated  in 
the  previous  two  sections,  either  through  use  of  the  tables 
in  Appendix  A,  through  direct  computation  using  Equation  5.1, 
or  from  the  charts  in  Appendix  B.   Further,  the  conditional 
probability  of  no  delay  given  no  initial  delay  is 


P[D  =  0|D.  =  0]   =    /   e  S   fv(x)dx 
4      L  0 


,    1/B    * 
"   (l/3  +  Xs)  ' 


The  equality  follows  from  the  assumption  of  Poisson  arrivals, 
with  rate  A  ,  of  serious  cases,  and  recognition  of  the 
integral  as  the  Laplace  Transform  of  the  Erlang  PDF  f(x). 
The  final  form  for  the  delay  probability  is  then 
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P[D>0]       =       1    -    (T-£ )        (1    -    P[T     >0]) 

i+  X  q 

as 


b.      Two   Server   System 

While    extension   to   allow  multiple   servers    is 
beyond  the   scope  of   this   thesis f    a  tractable    result    for   the 
case  where   there   are   only    two   servers    in   the   system   is 
presented  here,    although   the   early  part   of  the   derivation 
is    carried  out    for   general   c.      As   before, 

P  [D>  0]       =      1    -   P  [D  =  0]     . 

Conditioning  on  the  total  number  of  cases  in  the  system  at 
the  time  of  the  routine  case's  arrival  yields 


c-1 

P[D=0]       =         I      P[D=0|N=i]    P[N=i]  (5.2) 

i=0 


The  second  factor  is  merely  the  steady  state  system  size 
distribution  which  has  been  conveniently  tabulated  in  Ref. 
4  and  is  contained  in  Appendix  A.   Now  consider 


l 
P[D=0|N=i]       =         )"      P[D  =  0  |N     =  r,N  =  i]P  [N     =  r  |N  =  i].      (5.3) 

r=0 


In  this  expression,  the  second  factor  under  the  summation 
can  be  recognized  as  the  probability  of  r  routine  cases 
out  of  i  total  cases  currently  in  service.   Since  the 
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probability  that  any  one  case  is  routine  or  serious  is 
AR/(A  +Ag)  or  *s/(*ji+*s)  respectively, 


r  =   0  , 1 ,  .  . . ,  i ; 

i   =   0  , 1 ,  .  .  .  ,  c-1  * 


The  first  factor  is  more  difficult,  however,  and  requires 
further  conditioning  on  the  number  of  serious  case  arrivals 
during  a  routine  service  period: 


P[D=0|N    =  r,N=  i]       =  P[D=  0|A    =  j,N     =  r,N=  i]  (5.4 

R  j=0  S  R 


x    P[AS  =  j|NR=  r,N=  i]  . 


Since  arrivals  are  assumed  to  be  independent  of  the  status 
of  the  system, 


P[AS  =  j|NR=  r,N  =  i]   =   P[Ag  =  j]. 


Assuming  the  previous  distributional  assumptions  hold,  that 
is  serious  cases  arrive  in  a  Poisson  fashion  with  rate  \ 
and  the  Erlang  (k,8)  distribution  is  appropriate  to  describe 
the  service  times, 
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PtA  =j]   =    /   P[A  =j|Y  =  x]  fv(x)dx 


Since 

-X    x  (Ax)1' 
P  [A  =  j  |Y  =  x]   =   e 


S   Jl  3! 


and  fv(x)  was  given  in  Equation  (3.1)  , 


where 


P[As  =  j]   =   (k+;j  X)  pk  qj  (5.5) 


I                                                             x 
p   =   (y-2 )      and     q   =  (T-S^ ) 

J+AS  ^+As 


The  resultant  expression  can  be  recognized  as  a  Negative 
Binomial  PDF  although  that  fact  will  not  be  exploited  here. 
Since  this  derivation  is  restricted  to  the  case 
of  two  servers,  possible  values  of  i  in  formula  (5.4)  are 
limited  to  zero  and  one,  allowing  enumeration  of  all  possi- 
ble combinations.  For  i  =  0,  r  =  0  is  the  only  possibility 
Then 


P[D=  0|A^  =  j,N  =  0,N  =  0] 
S      K 


if   j  =  0,1 


(pk)j_1     if   j    2 
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because  if  no  serious  case  or  only  one  serious  case  arrives 
during  the  routine  case's  service  time,  his  service  will  be 
completed  uninterrupted.   If,  however,  more  than  one  serious 
case  arrives  during  his  service  time,  his  service  will  be 
uninterrupted  only  if  the  second  serious  case  arrives  after 
the  completion  of  service  of  the  first  serious  case  and  the 
third  serious  case  arrives  after  completion  of  service  of 
the  second  serious  case,  etc.   The  assumption  that  serious 
cases  arrive  in  a  Poisson  Process  at  rate  Aq  implies  that 
the  probability  of  a  serious  case's  service  completing  prior 
to  the  next  serious  arrival  is 

k         i    E    ^ 
p     =    (i ) 

since  a  service  time  is  the  sum  of  k  Exponential  phases, 
each  with  mean  duration  3.   Then  Equation  (5.4)  with  i  =  0 
and  r  =  0  becomes 

CO 

p[D=  0|N     =  0,N=  0]       =      P[A     =0]     +      [     (pk)  j~1    P[AQ  =  j]     (5.6) 
R  b  j  =  ]_  b 

Next  letting  i  =  1  and  r  =  0 


P[D=  0|Ag  =  j,NR=  0,N  =  1]       = 


if       j   =    0 


P[S  <  T]  (pk)  j    l  if       j 
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where  S  equals  the  time  to  completion  of  remaining  service 
of  the  active  serious  case  and  T  equals  the  time  to  the  first 
serious  case's  arrival.   This  expression  is  explained  as 
follows.   If  no  serious  case  arrives,  no  interruption  of  ser- 
vice occurs.   If,  however,  one  or  more  serious  cases  arrive, 
the  routine  case's  service  will  be  uninterrupted  only  if 
the  first  serious  case  arrives  after  completion  of  the  cur- 
rently served  serious  case's  remaining  service  time  and 
each  new  serious  case  arrives  after  the  completion  of 
service  on  the  previous  serious  case.   Therefore,  Equation 
(5.4)  with  i  =  1,  r  =  0  is 


P[D=  0|N     =  0,N=  1]       =      P[A     =0]    +    P[S<T]       £     (pk)  j-1P  [A_  =  j  ]  . 
R  S  jsel  b 


(5.7) 


Now,  letting  i  =  1  and  r  =  1, 


p[d=  o|as  =  j/NR=  i,n=  i; 


if  j  =  0 


(jP[S  >  T]+P[S  <_T]  (pk)  j  L  if  j  >  1 


Again,  no  service  interruptions  can  occur  if  no  serious 
case  arrives.   If  exactly  one  serious  case  arrives,  service 
can  continue  uninterrupted  only  if  the  other  routine  case 
either  remains  in  service  until  the  serious  case  arrives  and 
then  absorbs  the  resultant  preemption  (with  probability  one 
half) ,  or  if  he  completes  his  remaining  service  prior  to  the 
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serious  case's  arrival.   All  subsequent  serious  cases  must 
then  arrive  after  the  preceeding  serious  case  has  completed 
an  entire  service  period.   Therefore,  when  i  =  1,  r  =  1, 
Equation  (5.4)  becomes 


P[D  =  0 1  NR=  1,N  =  1]       =      P[Ag=0]     +    (jP[S>T]  (5.8) 


+    P[S  <  T])        I        (pk)D_1P[A0  =  j] 
j-1  S 


It    can  be    shown  that 


k,  j-1 


I       (p")JXP[A=j]       =      i— _-l  (5.9) 

j=l  S  (l-pkq)K 


and 


P[S^T]       =      ^-[1    -    p*]     .  (5.10) 


The  derivations  are  included  as  Appendix  C. 

The  means  now  exist  to  calculate  the  delay  proba- 
bility for  a  routine  case  in  a  two  server  system.   First, 
however,  it  will  be  convenient  to  adopt  the  following  simpli- 
fying notation: 


tt0      =      P[N=0] 


tt        =      p[n=  1] 
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Q0      =      P[A    =0J       =      pk 


Qi    -      I  (pk>j_1PtA  -j]     =    — 1         -  1 

3  =  1  °  (1  -  p    q) 

Q2      =      P[S<T]       =      ^-(1    -    pk) 


Then    for    c   =    2,    Equation    (5.2)    yields 


P[D=0]       =      P  [D=  0|N  =  0]7TQ    +    P  [D  =  0  |N  =  1]^ 


Now,    from  Equations    (5.3)    and    (5.6)  ,    it's    easy   to   see   that, 
for   the  two    server   case 


P[D=  0|N  =  0]       =      QQ    +   Ql 


and   further  that,    using  Equations    (5.3),     (5.7),    and    (5.8), 


P[D=0|N  =  1]       =        [Q0+QxQ2]q    +    {QQ    +     [j(l  -  Q2)   +Q2]Qx}p 


qQQ    +    (l-p)Q1Q2    +    pQQ    +    jpQ1    +    tPQ]_Q2 


QQ    +   |PQl    +    (i-lp)QlQ2 


yielding 


P[D=0]       =       (QQ  +Q1)tt0    +    (QQ    +  jpQ1    +    (1  -  jp)Q1Q2)T1     ( 
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Finally   then, 


P[D>0]       =      1    -    P  [D  =  0 


=   1  -  [(Q0+Q1)tt0  +  (QQ  +  ipQx  +  (1 -•|P)Q1Q2)tt1]  . 

(5.11) 


c.   Example 

The  following  example  is  provided  to  illustrate 
the  computational  procedure  required  to  obtain  this  delay 
probability.   Parameter  values,  though  hypothetical,  were 
chosen  so  as  to  provide  realistic  estimates  of  those  encoun- 
tered in  the  course  of  this  thesis. 

X    =   1.0/hr.       X    =   .0  5/hr.       8   =   .5  hrs . 
k   =   2      c   =   2      y   =   r-r-=l      p  *     — -  =   .5 

The  values  of  the  components  of  Equation  (5.11)  are  provided 
be  low . 

1 
p   =   j*-£ =   .9756 

e +  xs 

TTn   =   .3308   (from  Appendix  A) 


tt.   =   .3383   (from  Appendix  A) 
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QQ       =      pk      =       .9518 


CO 

Ql      =       .1    (P    ):"      PIA     <  j]       =       .0005 


j=l  S 


Q2      =      P[S<T]       =       .9640 


Direct    substitution   yields 


P[D>0]       =      1    -    {[(.9518    +    .0005)  ](  .3303)  +     [.9518    +   |(.9756) 


x    (.0005)    +    (1    -    j(. 9756)  )  (.0005)  (.9640)  ](  .3383)  } 


=       .3629 

It  is  interesting  to  compare  this  result  with 
that  obtained  using  the  same  parameters  but  considering  only 
a  FIFO  priority  discipline.   The  delay  probability  for  the 
single  priority  case,  as  taken  from  the  charts  in  Appendix 
B,  is  only  .1750,  substantially  less  than  that  for  a  routine 
case  in  a  two  priority  system.   But  this  is  to  be  expected, 
as  it  is  a  reflection  of  the  opportunity  for  additional  delay 
introduced  by  the  possibility  of  bumping  in  the  queue  and 
preemption  from  service. 
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VI.       CONCLUSIONS 

A.  DISCUSSION 

The  results  presented  in  this  thesis  provide  the 
machinery  to  conduct  a  sound  comparative  analysis  of  the 
merits  of  differing  boat  crew  allocation  schemes.   Investi- 
gation of  Tables  III,  IV,  V,  and  VI  shows  quite  dramatically 
the  improvement  in  the  various  measures  of  system  performance 
brought  about  by  increasing  the  number  of  available  servers. 
These  results  are  obtainable  without  great  computational 
capability  or  familiarity  with  the  underlying  mathematical 
theory.   The  use  of  one  or  the  other  of  the  two  methods  pre- 
sented is  a  matter  of  user  preference.   In  addition,  these 
techniques  suggest  the  development  of  system  standards  by 
which  proposed  crew  allocations  could  be  deemed  acceptable 
or  unacceptable.   Such  a  project  would  be  a  more  appropriate 
task  for  staff  members  of  the  Search  and  Rescue  Branch  at 
Coast  Guard  Headquarters. 

B.  RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

This  thesis  leaves  some  fairly  obvious  avenues  for  future 
work,  the  first  being  the  extension  to  the  general  case  of 
c  servers.   Perhaps,  by  first  allowing  three  servers,  the 
results  might  not  only  be  tractable  but  could  suggest  an 
extension  to  allow  any  number  of  boat  crews.   Another  area 
for  further  work  would  be  to  increase  the  number  of  priority 
classes.   As  previously  discussed,  boat  crews  are  actually 


51 


assigned  according  to  a  continuous  priority  scheme,  with 
each  arriving  case  being  assessed  as  to  its  severity  relative 
to  other  active  cases.   While  the  two  class  case  studied  in 
this  thesis  is  an  improvement  over  a  simple  FIFO  discipline, 
it  still  represents  a  gross  simplification  and  could  be 
improved  upon.   Finally,  the  computational  techniques  pre- 
sented herein  raise  the  possibility  of  the  creation  of  an 
interactive  computer  package  designed  to  supply  some  or  all 
MOE's  at  the  touch  of  a  computer  keyboard.   Though  not  com- 
plicated individually,  the  sheer  volume  of  computations  make 
hand  calculation  tedious  at  best. 

C.   PRIORITY  CONSIDERATIONS 

Considering  the  utility  of  approximating  the  actual 
priority  discipline  employed  by  an  RCC  controller,  and  of 
applying,  consequently,  the  methods  proposed  herein,  the 
point  of  view  of  the  SAR  planner  must  be  considered.   The 
consequence,  in  terms  of  loss  of  life  and  property,  of  a 
failure  on  the  part  of  the  Coast  Guard  to  provide  a  ready, 
available  rescue  response  in  the  event  of  a  serious  case 
is  far  more  grave  than  for  that  of  a  routine  case.   A  fishing 
boat  out  of  gas  in  a  protected  bay  can  tolerate  an  hour's 
delay  far  better  than  a  vessel  on  fire  or  taking  on  water 
and  in  danger  of  sinking.   Keeping  this  in  mind,  and  observing 
from  the  data  that  cases  of  a  serious  nature  are  the  rare 
exception  rather  than  the  rule,  the  two  priority  class 
assumption  leads  to  a  conservative  estimate  of  probable 
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system  performance  from  the  point  of  view  of  the  truely 
serious  case.   In  other  words,  the  actual  likelihood  that  a 
distress  case  requiring  rapid  emergency  response  would  not 
receive  immediate  attention  will  always  be  less  than  that 
predicted  by  the  queueing  techniques  described  in  this  work. 
This  conservative  posture,  then,  provides  an  upper  bound  or 
worst  case  prediction  for  the  serious  arrival,  and  therefore 
justifies  its  implementation  until  such  time  as  continuous 
priority  considerations  might  be  properly  integrated  into 
the  solution  methods. 
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APPENDIX  A 


TABLE  I 


STEADY-STATE  DISTRIBUTION  AND  MEAN  OF  NUMBER  IN  SYSTEM 
FOR  VARIOUS  M/E,/c  SYSTEMS 


Notation; 

N 

STATE  I 

P(N=  I) 

P(N  <  =1) 

RHO 

P( DELAY) 

L( GIVEN  K) 
LQ( GIVEN  K) 

(LQ  FOR  K  =  1) 
RATIO 


=   the  number  of  customers  in  the  system 
(including  those  being  served)  in  steady 
state  (a  random  variable) . 

=   the  condition  of  having  N  =  I. 

probability  that  N  =  I 

=   P(N  <.  I)  . 

=   p,  the  traffic  intensity. 

=   P(N>c)  =  P(T  >  0)  ,  the  probability  that 
~        q 
an  arriving  customer  cannot  begin  service 

immediately . 

=   L  =  E[N],  the  expected  number  of  customers 
in  the  system. 

=   L   =  E[N]  -  cp ,  the  expected  number  of 

customers  in  the  queue  (excluding  those 
being  served) . 

=   L  (k  =  1)  ,  the  value  of  L   for  the  corres- 

q  q 

ponding  M/M/c  system. 
=   L  /(L  (k=l)). 
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APPENDIX  B 
DELAY  PROBABILITY  (M/M/c  QUEUE) 

The  following  two  charts  provide  the  delay  probability 
P  [T  >  0]  for  the  M/M/c  queue  and  consequently  provides  an 
excellent  approximation  for  the  M/E,/c  queue  as  discussed 
in  Chapter  V. 
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APPENDIX    C 

1.      Derivation   of  P[S<_T] 

S      =      Excess  Life   of   customer   in   service 
T      =      time   to   arrival   of   next   serious    case 

The    PDF,    in   steady   state,    of   the   Excess    Life   of   a   non-negative 
random  variable  with   CDF   F   is    shown  by    Ross    [Ref.    13]    to  be 


e  mean 


Given  the  Erlang  distribution  as  the  distribution  of  service 
time 


fe(x)   =   j~-[l  -  F(x)] 


.  ,  ,l.i+l      1 

i  V  ¥ i   -jx 

=     k     l     —T< —  x  e 
Ki=0    1' 


1  k_1 
=   |  I   f(x) 

i=0   Erlang (i+1,3) 


Then 


1  k"X 

Vx)   =   £  I      F(x) 

i=l   Erlang(i+l,S) 


Now 


98 


-Ax 
P[S£T]       =         /      P[S<_T|T  =  x]     A      e      a      dx 
0  S 


Letting  T   =   x   and  using   the    independence   of  S    and  T   gives 


-Ax 
P[S<_T]       =        /      P[S<x]    A      e  dx 


00  -Ax 

/      F    (x)     Ac    e      b      dx 
0  e  S 


where 


*>    ,    k-1  -Ax 

/      £      [      F(x)  A      e      b      dx 

0        K    i=0      Erlang(i+1,B) 

k-1    ,         °°      °°         (3-x)  J      -s-x  -Acx 

I      b     /         I        -V-  e    3      AQ   e      S      dx 
i=0    K    0      j=i+l       J"  b 


1-0    D=i+1    (l  +  x    ,3+1    0J  J' 


i=0    j=i+l      i+As  i+As 


,    k-1 

i=0    j  =  i+l 
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r      j         i+1      r     j 

l        p  q  =  P    q  Z   P 
j=i+l  j=0 


i+1    1 
■  P    q  I=p 


i+1 

=  P 


Then 


[S<T]   =  i  I      p1+1 
i=0 

1   ,1  -  P  \ 


-  I  |  <!-pk> 


=  M(1  _  pk) 


2.   Derivation  of  I    (p^)3'1   P [A  =  j ] 

j-1  b 


v  /  k.  j-1  _r_    .,      r  ,  k.j-l,k+j-l.   k   i 
I    (p  ) J    P[A  = ]]   =    I    (p  ) J   (   ^   )  P  q 

j-1  s         j=l  : 


v  , k+ j -lx  ,    k  x  j 
I    (   h   )   (p  q)  J 

j-1    D 


Letting  9  =  p  q 


100 


,kxj-l 


i  (p*)^  p[A  =j]  =   i  (^-x) 


k+j"^     flJ 


j-l 


j-l 


[  I   (k+r1)9j(i-e)k-(i-e)k] 


(1-6)       j=0         : 


■i-TrU    -    d-6)k] 


(i-e) 


d-e) 


-  1 


(1-p    q) 


-    1 
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